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Abstract.  The  broadening  and  shift  parameters  for  a  number  of  Fe  I  lines  perturbed  by  atomic  hydrogen  are 
computed  using  the  interatomic  potential  due  to  Hindmarsh  etal.  (1967.  1970).  It  is  also  shown  that  the 
rms  radius  and  the  effective  radius  of  the  radiating  atom,  which  determine  the  force  constants  in  the 
interatomic  potential,  can  be  simply  related  each  other,  depending  on  the  orbital  quantum  number  of  the 
atomic  level. 


1.  Introduction 


The  collisions  with  neutral  hydrogen  in  the  solar  photosphere  broaden  and  shift  the 
absorption  lines.  While  the  pressure  shift  is  commonly  neglected  in  the  line  synthesis 
(Beckers  and  De  Vegvar,  1978),  the  pressure  broadening  compete  with  the  velocity 
broadening  (granulation,  waves,  and  turbulence)  in  determining  the  tvings  of  medium 
and  strong  lines.  The  classical  formula  used  to  compute  the  pressure  broadening  by 
neutral  particles  is  based  on  the  Van  der  Waals  attraction  between  the  colliding  atoms 
(e.g..  Unsold,  1955).  The  broadening  values  obtained  in  this  way  appeared  to  be 
generally  too  low,  about  a  factor  2,  in  comparison  with  the  experiment  and  with  the 
values  inferred  from  the  analysis  of  the  solar  lines  (Holweger,  1971a,  b). 

Hindmarsh  et  ai  (1967,  1970)  extended  the  classical  Van  der  Waals  broadening 
considering  both  attractive  and  repulsive  terms  in  the  interaction  potential  between 
radiating  and  perturbing  atoms.  O’Neill  and  Smith  (1980a,  b)  have  shown  that  this 
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semi-empirical  lheor\’  predicts  the  broadening  of  a  number  of  Ca  i  and  Fe  i  lines  by  noble 
gases  generally  with  an  underestimate  of  only  20“„. 

We  used  the  theory  of  Hindmarsh  et  al.  to  compute  the  values  of  the  damping  and 
shift  for  a  number  of  Fe  i  lines  produced  by  neutral  hydrogen  perturbers  at  temperatures 
typical  of  the  solar  photosphere.  Some  of  the  broadening  parameters  computed  in  this 
note  have  been  already  used  for  the  line  synthesis  in  the  presence  of  velocity  fields  (e.g.. 
Gomez  et  al.,  1987). 


2.  Method 


The  damping  constant,  or  half  width  at  half  intensity  maximum,  y,  and  the  shift,  /?,  in 
angular  frequency  unit,  appear  in  the  Lorenlzian  line  profile 


l{0i)  = 


y/?t 

(w-  Wq-  ' 


(1) 


where  oj  is  the  angular  frequency  of  the  radiation  and  Wq  refers  to  the  line  center. 

According  to  Hindmarsh  eta/.  (1970),  the  damping  •/,  and  the  shift />,  per  unit  number 
of  perturbers  are 


y,  =  ^  =  4r:  I  C»  !  5(0,  0 ,  (2) 

=£.  =  2tt  ^  1;o.6|qio.4  5(^^)  (3) 

N  \  8  / 


Here  v  is  the  relative  velocity  of  the  colliding  atoms  and  hC^  is  the  differential  force 
constant  of  the  term  r '  *  in  the  potential  of  interaction  of  two  atoms  separated  by 
distance  r.  The  factors  B  and  S  have  the  values  0.3  and  -  0.44,  respectively,  in  the  case 
of  a  pure  van  der  Waais  potential  (r~*);  in  general  they  are  the  integrals 


00 


5(0,0=  J  A:sin^[^(ojc  -  x'^)]dx , 

0 

(4) 

oo 

5(0,  ij)  =  J  X  sin(ax””  -  -  x"*)dx , 

(5) 

0 

and  the  parameters  o  and  ^  are 

0  =  0.5395'^  |C,;Jl/|Q|"•^ 

(6) 

^  =  0.780i;''''lCg|/|Q|■•^ 

(7) 

where  ACg  and  AC, 2  are  the  force  constants  of  the  terms  r"*  and  r~ 

respectively. 

in  the  interatomic  potential.  The  contribution  of  the  term  r  *  is  negligible  in  the  case 
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of  light  perturbers.  In  the  solar  photosphere  the  main  perturber  is  hydrogen,  then  we 
assume  ;  =  0. 

I  he  expression  of  Q  (Unsold,  1955)  is 


(8) 


where  a  is  the  dipole  polarizability  of  the  perturbing  atom,  which  has  the  exact  value 
of  0.67  X  10'  -■*  cm-\  and  r~  and  rf  are  the  quantum  mechanical  average  values  of  r- 
for  the  upper  and  lower  states  of  the  transition,  respectively.  Accurate  values  of  the  mean 
square  radii  for  all  assigned  terms  in  the  Fcl  energy  diagram  were  calculated  by  Warner 
(1969)  and  we  used  these  values  in  our  calculations.  Hindmarsh  etal.  (1967)  suggest 
for  C[2  the  expression 

AC,2  =  0.9  X  10 -'«(«, +  7^^)‘^  (9) 


TABLE  I 

Lines  and  atomic  parameters 


X 

(A) 

•‘I 

n; 

k 

£. 

W) 

(»2) 

!*«// 

(oo) 

5074.75 

1.923 

3.319 

1 

2 

23.77 

124.90 

22.59 

5137.39 

1.912 

3.226 

1 

2 

23.77 

145.20 

23.70 

5198.72 

1.548 

2.033 

0 

1 

11.36 

19.76 

10.38 

5217.40 

1.704 

2.426 

1 

0 

17.25 

94.38 

16.41 

5250.22 

1.323 

1.585 

0 

1 

9.83 

13.48 

9.24 

5250.65 

1.704 

2.426 

0 

1 

11.36 

19.76 

10.38 

5263.31 

1.714 

2.444 

1 

0 

17.25 

94.38 

16.41 

5302.31 

1.717 

2.444 

1 

0 

17.25 

94.38 

16.41 

5379.58 

1.799 

2.777 

0 

1 

9.73 

23.89 

11.02 

5395.22 

1.985 

3.432 

1 

2 

25.23 

140.90 

23.48 

5434.53 

1.405 

1.719 

0 

1 

11.78 

17.25 

9.95 

5506.79 

1.403 

1.709 

0 

1 

11.78 

17.25 

9.95 

5576.10 

1.745 

2.462 

1 

0 

18.09 

94.38 

16.41 

5635.83 

1.934 

3.077 

1 

1 

17.30 

196.40 

23.43 

6240.65 

1.548 

1.920 

0 

1 

11.36 

17.30 

9.96 

6246.33 

1.780 

2.426 

1 

0 

20.04 

94.38 

16.41 

6252.56 

1.574 

1.968 

0 

1 

9.77 

26.11 

11.34 

6297.80 

1.548 

1.916 

0 

1 

11.36 

23.01 

10.89 

6301.51 

1.790 

2.444 

1 

0 

20.04 

94.38 

16.41 

6302.50 

1.797 

2.462 

1 

0 

20.04 

94.38 

16.41 
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where  the  numerical  constant  (in  erg)  was  fixed  empirically  and  and  are  the 
distances,  from  the  nucleus  of  the  radiating  and  perturbing  atoms,  respectively,  at  which 
the  unperturbed  radial  charge  density  has  the  value  0.012  atomic  units.  We  used  the 
Coulomb  approximation  (e.g.,  Seaton,  1958)  to  compute  the  effective  radius,  since  its 
definition  involves  the  knowledgement  of  the  atomic  wavefunction  at  long  distances 
from  the  nucleus. 


3.  Results 

In  Table  1  we  list  the  20  Fe  i  lines  whose  broadening  and  shift  parameters  were  computed 
for  this  note.  For  each  line  we  report  the  effective  principal  quantum  number,  n*,  and 
the  orbital  quantum  number,  /,  of  both  energy  levels,  which  are  needed  to  specify  the 
corresponding  Coulomb  wavefunction.  The  6th  and  7th  columns  of  the  table  contain 
the  mean  square  radius  of  the  lower  and  upper  level  respectively,  derived  from  Warner 
(1969).  The  last  column  contains  the  values  of  the  effective  radius  of  the  radiating  atom, 
R^ff.  The  effective  radius  for  atomic  hydrogen  is  4.38  atomic  units. 

In  Table  II  we  list  for  each  line  the  multiplet  number,  the  level  notation,  the  damping 


TABLE  II 

Broadening  and  shift  parameters 


A 

(A) 

Multiplet 

No. 

level 

lowtr 

level 

upper 

l/N 

(10“' rod  cm®s~‘) 

Tv.twIN 

(10“*  rodem'a"*) 

0IN 

(m/a) 

5074.75 

1094 

if’F* 

e’Gs 

6.06 

2.30 

64.0 

5137.39 

1090 

A' £>4 

6.59 

2.48 

70.9 

5198.72 

66 

y'f’a 

1.59 

0.85 

14.8 

5217.40 

553 

«'i?3 

3.31 

2.07 

26.8 

5250.22 

1 

o' Do 

z'’D, 

1.36 

0.61 

13.7 

5250.65 

66 

o'Pj 

y'Fa 

1.59 

0.85 

15.0 

5263.31 

553 

a'Xlj 

e'K, 

3.31 

2.07 

27.0 

5302.31 

553 

a'£>, 

e'Dj 

3.31 

2.07 

27.2 

5379.58 

928 

a* 

1.73 

1.05 

14.5 

5395.22 

1143 

6.49 

2.43 

74.3 

5434.53 

15 

o'f, 

a®Z>o 

1.50 

0.72 

15.7 

5506.79 

15 

o'Fj 

a' £>3 

1.50 

0.72 

15.9 

5576.10 

686 

a'Fi 

e®£)o 

3.32 

2.06 

28.4 

5635.83 

1088 

x'Fi 

tt'F, 

2.58 

2.95 

24.3 

6240.65 

64 

o'P, 

a' Fa 

1.50 

0.74 

17.6 

6246.33 

816 

x'Fj 

e'£?3 

3.32 

2.04 

33.1 

6252.56 

169 

a=>ffa 

a'Gs 

1.80 

1.11 

17.4 

6297.80 

62 

o'Fi 

y'Oa 

1.71 

0.97 

17.8 

6301.51 

816 

x'Fa 

e'O, 

3.32 

2.04 

33.4 

6302.50 

816 

x'Ft 

e'£>o 

3.32 

2.04 

33.4 
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constant  (full  width  at  half  intensity  maximum)  and  the  shift  in  velocity  unit.  The 
damping  constant  is  calculated  at  the  temperature  T  =  5000  K  and  per  unit  number  of 
neutral  hydrogen  atoms.  To  change  the  temperature  it  is  sufficient  to  multiply  by 
(7/5000)“  '*  (O’Neill  and  Smith,  1980a).  As  a  matter  of  comparison,  we  report  also  the 
broadening  computed  with  the  pure  Van  der  Waals  theory.  It  appears  that  the  Van  der 
Waals  values  are  generally  smaller  than  the  values  calculated  by  including  the  term  r  ' 
in  the  interatomic  potential,  which  confirms  the  results  of  O’Neill  and  Smith  (1980b). 
A  comparison  with  the  euipirical  damping  constants  determined  by  Simmons  and 
Blackwell  (1982)  for  lines  5434.53,  6252.56  and  multiplet  62  indicates  that  our  data  are 
slightly  larger  ( ~  25  % )  than  the  data  based  on  the  fit  of  solar  lines.  The  shift  parameter 
has  been  computed  at  the  temperature  of  5000  K  and  for  a  perturber  number  density 
of  10'^  neutral  H  per  cm^.  This  density  corresponds  to  a  depth  between  0  and  50  km 
in  the  VALC  atmospheric  solar  model.  Since  the  hydrogen  number  density  decreases 
exponentially  with  depth  in  the  photosphere,  the  values  we  found  for  the  pressure  shift 
of  Fei  lines  confirm  the  findings  of  Beckers  and  De  Vegvar  (1978)  that  for  the  Fei  lines 
the  pressure  shift  is  negligible  with  respect  to  the  convective  blueshift  produced  by  the 
solar  granulation. 
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Fig.  1.  Ratio  of  the  effective  radius,  R^,  to  the  r.m.s.  radius,  of  an  atomic  level  versus  the  r.m.s. 
radius.  Circles  refer  to  / »  0,  crosses  to  /  -  1 ,  and  squares  to  /  -  2,  with  /  being  the  orbital  quantum  number 
of  the  level.  Both  radii  were  computed  on  the  base  of  the  Coulomb  approximation  to  the  atomic  wave- 

function. 
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Finally.  Figure  1  shows  the  relation  e.xisting  between  the  effective  radius  and  the  mean 
square  radius  of  a  level  when  both  are  computed  on  the  base  of  the  Coulomb 
approximation.  Usually  the  ratio  of  the  effective  radius  to  the  r.m.s.  radius  is  approxi¬ 
mated  by  2  (e.g..  Hart,  1974).  However,  Figure  I  shows  that  this  ratio  is  a  raonotonic 
decreasing  function  of  the  r.m.s.  radius  which  depends  on  the  orbital  quantum  number 
of  the  level.  This  relation  provides  a  quick  way  to  estime  the  effective  radius,  e.g.,  by 
using  the  Unsdld  (1955)  formula  to  evaluate  the  mean  square  radius. 
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